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Abstract 

Top-quark physics plays an important role at hadron colliders such as 
the Tevatron collider at Fermilab or the upcoming Large Hadron Col- 
lider (LHC) at CERN. Given the planned experimental precision, de- 
tailed theoretical predictions are mandatory. In this article we present 
analytic results for the complete electroweak corrections to gluon in- 
duced top-quark pair production, completing our earlier results for 
the quark-induced reaction. As an application we discuss top-quark 
pair production at Tevatron and at LHC. In particular we show that, 
although small for inclusive quantities, weak corrections can be size- 
able for differential distribution. 



Heisenberg Fellow 



I. Introduction 



Top-quark physics plays an important role at the Tevatron and will be an equally im- 
portant topic at the upcoming LHC. In view of the large production rate, amounting 
to (10^) top-quark pairs for an integrated luminosity of 200 fb^\ precise and direct 
measurements will be possible, which require a similarly detailed theoretical under- 
standing of these reactions. Both single top-quark production as well as top-quark pair 
production have been studied extensively in the past. The differential cross section 
for top-quark pair production is known to next-to-leading order (NLO) accuracy in 
quantum chromodynamics (QCD) dElElllllSll. In addition, the resummation of log- 
arithmic enhanced contributions has been studied in detail in Refs. L6..7. .8..9. irOl[TT1l . 
Recently also the spin correlations between top-quark and antitop-quark were calcu- 
lated at NLO accuracy in QCD Ed. 

Although formally suppressed through the small coupling, weak corrections can also 
be significant due to the presence of large Sudakov logarithms (see e.g. Refs. 11131 [T?l 
and references therein), which were also studied in the context of y- and Z-production 
at hadron colliders llfSl fT6l [TTll . The origin of these large logarithms is easily un- 
derstood: at high scale the massive gauge bosons W and Z behave essentially like 
massless bosons. CoUinear and soft phenomena thus lead to large negative correc- 
tions. In strictly massless theories like QED or QCD these contributions are cancelled 
through similar positive terms from the real corrections. This cancellation does not 
take place in the case of weak interaction because the real and virtual emission leads 
to different experimental signatures. Given that top-quark pair production at high scale 
is an ideal tool to search for new physics it is clear that the precise knowledge of the 
weak corrections in this region is of paramount importance. In Ref. lITHll electroweak 
corrections to top-quark pair production in hadronic collisions were investigated for 
the first time. More precisely, the partonic sub-processes qq tt and gg tt were 
studied. In a subsequent study [19| parity violating asymmetries were analysed in a 
two Higgs doublet model and the minimal supersymmetric standard model. 

In the original works some contributions were omitted. For the quark-antiquark ini- 
tiated process the gluon-Z box contributions (Fig. |lT] a) and the corresponding real 
corrections (Fig. |LT]b) are missing in Ref. |18|. They were recently evaluated in 
Refs. II20I I2TI . For the gluon fusion process a class of contributions related to tri- 
angle diagrams (Fig.|rT]c) are missing in Ref. [T8|, as noted also in Ref. [22J, where 
the calculation of Ref. [ 1 8 1 has been repeated for the gluon induced top-quark pair 
production. In Ref. |22| no analytic results are presented. In view of the importance 
of the analysis for upcoming experiments it is the purpose of this paper to repeat the 
original calculation ifTHI — including all missing contributions — and present compact 
analytic results well suited for the experimental analysis. 

The outline of the paper is as follows: In section HII we present the calculation of the 



1 




virtual electroweak corrections to top-quark pair production through gluon fusion. In 
contrast to the contributions from quark-antiquark annihilation there are no real cor- 
rections contributing to the weak corrections. The virtual corrections are thus infrared 
finite and represent the complete weak corrections for this channel. Compact analytic 
expressions in terms of scalar one-loop integrals are given in section HII and in the ap- 
pendix. In section Uni we present numerical results for the gluon fusion process at the 
parton level. Furthermore we combine the gluon channel with the quark-antiquark 
annihilation process (with the elctroweak corrections taken from Ref. L2QJ), fold them 
with parton distributions and give results for the corrections to the total cross section 
and to pt- and M^f- distributions relevant for the Tevatron and the LHC. 



II. Electroweak corrections to gluon fusion 
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Figure II. 1: Bom diagrams for top-quark pair production via gluon fusion. 



To setup our notation we start with the QCD tree-level contribution. The three contri- 
butions to the amplitude are shown in Fig. III. 11 Evaluating the Feynman diagrams we 
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obtain the well-known leading order differential cross section 



dz A^(l-|32z: 



2,2^2 



(l - PV + 2P2(1 - p2)(l (n.l) 



where is the number of colours, a, the strong coupling constant and P the velocity 
of the top-quark in the partonic centre-of-mass system: 
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(s denotes the partonic centre-of-mass energy squared). The cosine of the scattering 
angle is denoted by z. Here and in what follows it is convenient to use the abbreviation 
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A factor 1/(4(^2 _ i)2^ ^^.^^ 

averaging over the incoming spins and colour is included 

in the result above. 
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Figure II. 2: Sample diagrams for the virtual corrections. F stands for all contributions 
from gauge boson, goldstone boson and Higgs exchange. 



For the calculation of the next-to-leading order weak corrections we use the 't Hooft- 
Feynman gauge (i?^-gauge) with the gauge parameters ^' set to 1. The longitudinal 



3 



= ^ dn ^. - - . , ■ (11.4) 



degrees of freedom of the massive gauge bosons Z and W are thus represented by the 
goldstone fields % and (|). Ghost fields do not contribute at the order under consider- 
ation. Sample diagrams are shown in Fig. III. 21 The Cabibbo-Kobayashi-Maskawa 
mixing matrix is set to 1 . 

Before presenting the results for the weak corrections, let us add a few technical re- 
marks. We use the Passarino -Veltman reduction scheme ll23l to reduce analytically 
the tensor integrals to scalar integrals. For these the following convention is used: 

For the UV-divergent integrals we define the finite part for the one-point integrals Aq 
and the two-point integrals Bq through 

Ao(m^) = m^A-l-Ao(m^), 
5o(/,mim^) = A + 5o(/,m?,m^), (II.5) 

with A = 1/e — Y + ln(47r). The renormalisation is performed in the counterterm for- 
malism, where the bare Lagrangian L is rewritten in terms of renormalised fields and 
couplings: 

1/2 1/2 

= L{^'R,AR,mR,gR) + Lct{^R,AR,mR,gR). (n.6) 

The contribution L{^R^AR,mR,gR) gives just the ordinary Feynman rules, but with 
the bare couplings replaced by the renormalised ones. The complete list of Feynman 
rules can be found for example in Ref. |24|. The contribution Lct{^R,AR,mR,gR) in 
Eq. (III. 61 ) yields the counterterms, which render the calculation ultraviolet (UV)-finite. 
Some of the resulting diagrams are shown in Fig. III. 31 For the present calculation only 
wave function and mass renormalisation are needed. No coupling constant renormali- 
sation has to be performed. This is a consequence of the fact, that although the weak 
corrections appear as a loop-correction, they are still leading-order in the electroweak 
couplings. Mass and wave function renormalisation are performed in the on-shell 
scheme: 

"^/,o = mt + dnit, (11.7) 
= (^zf'^y^Vf'^= (^l + ^5zf'^^^f'^. (II.8) 

The renormalisation constants are thus given in terms of self-energy corrections E and 
their derivatives: 

5Zv = ^(5zf + 5zf) = -Ev(;?2^m,2)-2m,2^(Ev + Es) 

5Za = ^(5zf-5zf) = -EA(p2^m,2), 

^ = -Ev(/ = m,2)-Es(/ = m,2). (II.9) 
mt 
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Figure II. 3: Sample counterterm diagrams. 



The functions Ev,a.s can be found for example in Ref. ll24l . Here we need only 5Zv 
and dnit/nit. Their explicit form in terms of scalar integrals — using the notation of 
this paper — reads: 
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where sw (cw) denotes the sine (cosine) of the weak mixing angle. The vector {g[) 
and axial-vector (g'^) couplings of the top-quark to the Z-boson are given in terms of 
the weak isospin T.^ and the electric charge Qf for a fermion of flavour /: 
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The coupling of the top-quark to the W-boson is given by 

1 
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(11.12) 
(11.13) 

As usual a stands for the fine structure constant, which will be taken as running cou- 
pling evaluated at the scale 2mt. The mass of particle / is denoted by m, . The abbrevi- 
ations Bq'^'^ for the two-point scalar loop-integrals are defined in the appendix. Note 
that the photonic corrections form a gauge independent subset and are not included in 
Eqs. III. lOllII.l H and in the following discussion. 

For the following discussion it is convenient to separate the weak corrections into 
the contribution from vertex-diagrams {t-, w-channel, 5-channel, Fig. III. 21 a,d), self- 
energy-diagrams (Fig. III.2l b) and box-diagrams (Fig. III.2l c). The triangle diagrams 
in (Fig. III. 21 e) are finite without renormalisation and will be studied separately. The 
differential cross section at next-to-leading order is decomposed as follows: 
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We start with the analytic results for the triangle diagrams. For the Higgs and (Z + %)- 
terms we obtain 



do 



1 



dz 



mt^{s - 4m,^)Co + mb^{s - Amb^)Cl - 2{mt^ + nit' 



do 



A 



a 



-J— = 16-O0g«^^2(i_ 132^2) 
2 l^2 



/ 2/^? I b 2 ✓-'I 



' 



(11.16) 



(11.17) 



A factor 1/{4{N^ - 1)^) from averaging over the incoming spins and colour is again 
included. The integrals Cq' are defined in the appendix. As a consequence of Furry's 
theorem only the axial-vector induced terms contribute in the case where the Z-boson 
appears in the 5-channel. Furthermore, the Landau- Yang theorem forbids the decay of 
an on-shell vector boson into two identical massless on-shell spin-one bosons. There- 
fore, the poles from the propagators of the Z-boson and the % are cancelled in the 
02_f^-term, as evident from Eq. (III. 171) . 

For the remaining vertex corrections with a gluon in the i'-channel we obtain 
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The remaining contributions to the differential cross section are listed in appendix IHl 

Before showing concrete results for hadron colliders we first discuss several checks of 
our result. We performed two independent calculations yielding also two independent 
numerical computer codes. We checked that we obtain the correct structure for the 
UV singularities yielding a finite result after renormalisation. In our notation this 
corresponds to a stringent test of the coefficients of the Aq- and J5o-integrals. The 
behaviour of the Higgs corrections close to threshold and for light Higgs bosons is 
well understood (see for example Ref. [25 ( and references therein). This allows to test 
the Higgs contributions for a very light Higgs near threshold. For a ff-system produced 
through a vector current the Higgs correction is given by the factor (1 +//thr('")) with 
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Figure 11.4: Comparison of the full corrections from Higgs boson exchange with the 
corrections based on Eq. (III.23I) (mn = 0.1 MeV). 



T] defined through 

>l = £;i-i- (11.25) 

The line shows the Higgs contribution according to Eq. (III.23I) . The crosses show the 
result from the full Higgs boson induced correction. The enhancement proportional 
to l/P is well recovered by the calculation. The numerical agreement is between four 
(ri = 10^^) and two digits (Tj = 10^^). As we shall see later the enhancement close 
to threshold for a light Higgs can still be observed even for a Higgs mass of 120 GeV 
although reduced to 2%. A similar test has been performed for a light Z-boson. Again 
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we find perfect agreement. In addition we compared analytically the results for the 
t- and M-channel vertices, self-energies and Higgs triangle diagrams with those from 
Ref. [ 18|. We find complete agreement, after the correction of some typos in Ref. HTSll . 
Using the same input parameters we also compared the plots shown in Ref. [iTSi and 
found agreement. Finally we compared numerically with the results of Ref. ll26ll and 
found perfect agreement. However, we are in disagreement with the results published 
recently in Ref. lE^ . In particular for the pr and M^f distribution we find negative 
corrections close to threshold while the corrections are positive in Ref. Ii22il (Fig. 1). 
Our findings are also confirmed in Ref. L26J . 



III. Numerical results 



^ 5 



c 
o 







o 

(D 
i_ 
i_ 

o 
o 

(D 

CD 
> 

CD 



-5 
■10 
■15 



— 1 1 1 1 1 r- 



— 1 1 1 1 1 r- 




350 400 450 500 



10 



10 



Figure III.l: Different contributions to the electroweak corrections: Vertices (long- 
dashed), self-energies (dashed), boxes (dotted), triangles (dash-dotted). The sum is 
shown as full line. 



In this section we present numerical results for the gluon fusion process at order a^a. 
We use the following coupling constants: 

= -p^' a, = 0.1, sw^ = 0.231, 
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and, if not stated otherwise, the following masses: 

mz = 91.1876 GeV, miv = 80.425 GeV, m// = 120 GeV, 

nib = 4.82 GeV, nit = 172.7 GeV. 

The parameters sw, mz and mw are, of course, interdependent. Nevertheless, we expect 
that, using the MS value for the weak mixing angle, some of (uncalculated) higher- 
order corrections are included and, therefore, a better phenomenological description is 
achieved. The difference between using the MS and the on-shell input respectively, is 
formally of higher- order. 

In Fig. lIII.ll the separate corrections as defined in Eq. (III. 151) are shown at the parton 
level. We normalise the result to the leading-order gg — > tt process. The sum of the 
different contributions is shown as a solid line in Fig. IIII.ll For moderate partonic 
energies the corrections are of the order of a few percent as one might have expected 
for a weak correction. Near threshold the cross section is dominated by the triangle- 
and the box-diagrams. Both are of the same order of magnitude, but opposite in sign, 
leading to a significant cancellation. It is worth noting at this point that the inclusion of 
the (Z-|-%)-triangle diagrams — neglected in Ref. lITSl — decreases the result by about 
2% in the threshold region (m// = 120 GeV). The (Z-l-x)-term dominates the triangle 
contributions. In Fig. IIII.2l we illustrate the effect of <^z+i by comparing the full result 

with the result where o^^^ is neglected. Close to the threshold the aforementioned 
2% difference is observed. The (Z + %) -contribution accidentally compensates the 
positive contribution from Higgs exchange, which however, becomes small about 20 
GeV above threshold. For energies above 600 GeV the contribution from the (Z -t- 
X)-triangle becomes negligible. The behaviour of the box contribution close to the 
threshold, shown in Fig. IIII.ll is a consequence of the Higgs effects discussed before. 

On the other hand for a partonic centre-of-mass energy of around 500 GeV the large 
negative Sudakov logarithms start to become important and amount to more than 10 
percent at a few TeV If we compare the relative size of the weak correction for gluon 
and quark-antiquark induced reactions at large energies, we find that they are twice 
as large for the quark-antiquark process. This can be qualitatively understood by just 
counting the external lines which can emit W- and Z-bosons, and observing that the 
corrections start to be dominated by Sudakov logarithms. 

The dependence on the Higgs mass, i.e. the relative corrections for different Higgs 
masses are shown in Fig. lIII.3l For Higgs masses larger than 2m; we include the width 
of the Higgs boson in the ^-channel propagator. 



The corrections are strongly dependent on m// with a variation of nearly 6% in the 
threshold region. Let us now address the effects of the weak corrections on hadronic 
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Figure ni.2: Comparison between the electroweak corrections: with the Z and % trian- 
gle diagrams included (full) and without these contributions (dashed). 
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Figure III. 3: The dependence of the partonic cross section from the Higgs mass: nin = 
120 GeV (full line), mh = 180 GeV (dashed), rriH = 240 GeV (dashed-dotted) and 
niH = 1000 GeV (dotted). 
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Figure in.4: Electroweak corrections to top-quark pair production at the Tevatron 

(upper figure) and the LHC (lower figure) for three different Higgs masses {mn = 
120 GeV (full Une), mn = 200 GeV (dashed), mn = 1000 GeV (dashed-dotted)). 
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observables. We used the parton distribution function CTEQ6L [27| evaluated at a 
factorisation scale = Imj. With the input mentioned before we obtain 

^Tev _ og^ + o2^ = 4.18pb + 0.19pb = 4.37pb 

and 

o^HC = o^?c+o^«c = 56pb + 366pb = 422pb 

for the leading-order cross section at the Tevatron respectively LHC. The leading-order 
estimate is significantly smaller than the QCD corrected (NLO -i- resummation) result 
of about 6.7 pb for the Tevatron 1 1 1 1 and about 794 pb for the LHC [28 1. To some ex- 
tend the large QCD corrections are of universal character and it is plausible that lowest 
order and electroweak corrected (total and differential) cross sections will be affected 
by similar corrections. Therefore we will, in the following, only present their relative 
size. In a first step we study the weak corrections to the total cross section at the Teva- 
tron and the LHC as a function of for three different Higgs masses ( Fig. lIII.4l) . The 
results include both quark-antiquark annihilation and gluon fusion. As expected the 
corrections to the inclusive cross section amount to a few percent only. Most of the 
top-quark pairs are produced close to threshold, where the weak corrections at parton 
level amount to a few percent only. The relative corrections are essentially given by 
the threshold behaviour of the quark-antiquark channel for the Tevatron and the gluon 
channel for the LHC. The integrated cross section samples a wide range in Mfj, and 
the marked Higgs mass dependence of the partonic cross section close to threshold is 
washed out. The correction is nearly independent of the top-quark mass for values of 
m-t between 165 and 180 GeV. Given the experimental precision at the Tevatron and 
the LHC it is unlikely that the weak corrections can be seen in the total cross section. 
Differential distributions in pj and Mft are affected however, strongly. Indeed, the 
corrections to the Mff-distributions can be directly read off from Fig. IIII.31 as far as 
the gluon induced channel is concerned. Note that for the quark-antiquark process the 
situation is more involved due to the presence of real corrections ll20l . In principle 
it might be possible to establish the enhancement induced by the top-quark Yukawa 
coupling for a relatively light Higgs boson. However, the difference of 6% between 
a light (120 GeV) and a heavy (1000 GeV) Higgs boson could be masked by QCD 
uncertainties which are particularly large in the threshold region. In contrast, the weak 
corrections amount to more than ten percent at high energy when the Sudakov sup- 
pression becomes large. Therefore we study their effect on differential distributions 
at large momentum transfer. With the differential distribution in the top-antitop in- 
variant mass (M/f = a/ {h +kt)^) being a sensitive tool in the search for new physics, 
this question is of particular importance. At large momentum transfer two competing 
effects must be considered: The logarithmically increasing Sudakov logarithms, and 
the increasing statistical uncertainty. 

To get a rough idea about the possible statistical sensitivity we show in Fig. ini.5l and 
Fig. lIII.6l the leading-order differential cross sections in M^f respectively pr- For LHC, 
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Figure ni.5: Leading-order differential cross section for LHC as a function of pj and 
Mtf. Shown is the sum (full) and the contributions from gluon fusion (dashed) and 
quark-antiquark annihilation (dotted). 
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Figure ni.6: Leading-order differential cross section for Tevatron as a function of pj 
and Mtj. Shown is the sum (full) and the contributions from gluon fusion (dashed) and 
quark-antiquark annihilation (dotted). 
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about 10'' events are expected for an integrated luminosity of 200 fb^ and large values 
of Mjt and pj will be accessible. An interesting behaviour of the relative importance 
of gluon- versus quark-induced processes is observed for the LHC. For low Mtj and pr 
the gluon channel dominates. However for pj larger than 1 TeV the quark-antiquark 
annihilation process takes over, as a consequence of the change in relative importance 
of quark-antiquark versus gluon luminosities. For the Tevatron only moderate values 
of 300 GeV for pj and 700 GeV for M;^^ are accessible at best. Furthermore, the 
quark-antiquark annihilation process dominates through out. 

The relative corrections of the differential distributions in pj and M;^^ are shown in 
Fig. 1111.71 and Illl.Sl for the LHC and Tevatron. For large values of pj and M^^, ac- 
cessible at the LHC, sizeable negative corrections are predicted, reaching ten up to 
fifteen percent. In contrast, the corrections are far smaller at the Tevatron. Taking 
niH = 120 GeV, they are +3% close to threshold and -5% for M„- around 800 GeV, 
leading to a distortion of the differential distribution by 8%. It remains to be seen, 
whether QCD and experimental uncertainties can be pushed below this level. A rough 
guess of the statistical uncertainty expected for LHC and Tevatron can be deduced 
from Figs. IIII.9I and Fig. IIII.IOI The estimated number of events with Mfj > M^"', 
based on a luminosity of 200fb^ for the LHC (8fb^ for the Tevatron) is used to es- 
timate the statistical uncertainty and compared with the relative corrections, evaluated 
for the same sample. It will be difficult to observe the effect of the electroweak correc- 
tions at the Tevatron. At the LHC, with the large sample of top quarks, the statistical 
precision will match the size of the electroweak Sudakov logarithms, and eventually 
of the Higgs enhancement in the threshold region. 

Before closing this discussion, let us mention that also the dependence of the cor- 
rections on the bottom-mass was investigated and the full dependence on the bottom 
quark mass was kept throughout. Furthermore the results can also be used to study 
weak corrections for Z^-quark pair production. This topic will be discussed in a sub- 
sequent publication. For a massless bottom-quark the W, (|) contributions and the Z, 
X triangles differ at most one percent from the massive case. Hence the effect of the 
bottom-mass is negligible at hadron level. 



IV. Conclusion 

In this article the complete electroweak corrections to gluon induced top-quark pair 
production are calculated. In contrast to earlier publications all contributions of one- 
loop-order are taken into account and presented in the form of compact analytic ex- 
pressions — well suited to be used in the experimental analysis. Furthermore the full 
dependence on the bottom-quark mass is kept. This allows to calculate also weak cor- 
rections for bottom quark production using the results presented here. We have shown 
that the corrections are negligible for the total cross section. For differential observ- 
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Figure ni.7: The relative corrections to the pr- and M^f- distribution for the LHC for 
ruH = 120 GeV (bold line) and mn = 1000 GeV (thin line). 
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Figure ni.8: The relative corrections to the pr- and Mtf- distribution for the Tevatron 
for = 120 GeV (bold line) and m// = 1000 GeV (thin line). 
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Figure ni.9: The relative corrections to the pr- and M,f-distributions, integrated from 
a lower value in pr respectively M^f to the kinematic Umit, for the LHC and two 
Higgs masses (m^ =120 GeV (bold line), m// = 1000 GeV (thin line)). The smooth 
curve gives an estimate of the corresponding statistical uncertainty for an integrated 
luminosity of 200 fb-^ 21 





Figure in. 10: The relative corrections to the pr- and M,f-distributions, integrated from 
a lower value in pj respectively M,f to the kinematic limit, for the Tevatron and two 
Higgs masses (m// =120 GeV (bold Une), niH — 1000 GeV (thin line)). The smooth 
curve gives an estimate of the corresponding statistical uncertainty for an integrated 
luminosity of 8 fb~^ . 22 



ables like the -distribution or the distribution in the invariant tt-mass the corrections 
can be sizeable. In particular we find corrections up to fifteen percent in kinematic 
regions which are accessible at the LHC. 
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Dl{z) 
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(A.l) 



To evaluate the Higgs triangle contribution for mn > 2mt we give the result for the 
corresponding vertex integral: 



Co(0,0,^,m2,m2,m2) = j- 



1 
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l+P. 
1-|3„ 



(A.2) 



with 



Arn^ 



(A.3) 



B. Analytical results 



The corrections not yet listed in section|niare divided in self-energy (Fig. III. 21 b). vertex 
(Fig. III.2l a) and box (Fig. III.2l c) corrections. Self-energy corrections: 
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+ (gf ( 1 + 7p2 - p4 - 2p6 + (3 + 8p2 - 4p4)pz 
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(1 + - p4)p V + P^z^) - (3 - + 5P^ - 

(3 - 8p2 + I2p4)[3^ + 3(1 _ [32 _ (34^(32^2 _ pS^S)^ j^z^^) 

(z^-z), 



(B.ll) 



+ 



4p(l-z2)(4p-z-2pz2) . 



5(l + p2 + 2pz) 



(^Ao(mi,2)-Ao(mM.2)) 



1 

2^ 



4{mw^ - mb^)z{5 - - (3 - 2^^)z^) 



- 5 ( 4p - (5 + - 4p^)z - 4pz^ + (3 + - 2p4)^3 



^0 



+ 



1 



4{mw^ — rrih^) 



25p(l + pz 

(6p3 + (5 - 8p2 +4p4)2+ (5 - 12p2)p22 _ (3 _4[32 + 2^4)^3 

(3 -4p2)p^4^ +5(2(5 - p2)p3 + (1 _ p2)((5 ^ I2p2 _4p4)^ 



+ 9pz^-(3 + 4p^-2pV)-P(3 + 5p")z' 



i2\4 



w 



{4(mw^ - mb^)f + 8mv7^5(2 + P^ + Pz) - 8mi^5P(2P + z) 



1 

+ 2s 

+ 52(7 + 2p2 + pV2(l + p2)pz + 2pV) 
1 



+ 



85P 



- {A{mw^ - mb^)fz{5 - %f - (3 - 2p2)z2) 
+ 8mw^5(2(l+p2)p-(5-5p2-8p4)z 

- 2(l-p2)Pz2+(3 + 3p2-2pV) 

- 8mfc2^(2(l+p2)p-(l + p2)(5-8p2)z 

- 2(l-p2)Pz2 + (l + p2)(3-2pV) 

- / ( _ 4( 1 + 4p2 ^ p4)p ^ (5 _ 3Qp2 ^ p4 _ 8p6)^ 

+ 4(1 + 2p2 - p4)pz2 - (3 +4p2 + 1 lp4 - 2p6)z3)] 

1 r 

^(l+pz)(l + p2 + 2pz) . 

4(miy2 _ ;„^2)p(^ + p) ( 1 _ 3p2 ^ ( 1 _ 2p2) p^ ^ 2p V + p323) 

5( 1 + p2 + 2pz) (2 + 5p2 - p4 + (3 - 6p2 + 2p4) pz 



+ 
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(7 - 2p2)pV + (2 - p2)p3^3 _ ^\A]b^(,) 



1 

- 8mfc2^P(P + z)(2 + Pz) 

+ ^2(1 + p^) + 2p2 + p4 + 2(1 + p2)Pz + 2pV)] C^(z) 

lOmiv"^ + 2m^4(l + 2pV) - 4miy V^(3 + P^z^)) 
+ 4*2 (^miv2^ J 1 + 8p2 + 3p4 + 4(3 + 2p2)pz + 6p V) 

- m/,2(ll-4p2 + 3pV8(l + p2)Pz + 2(6 + p2)p2z2)^ 

+ ^3 ( 1 + p2 + 2pz) (7 + 2p2 + p4 + 2( 1 + p2) pz + 2p V)] < (z) 
+ (z^-z), 



1 

+ is 



+ 



(B.12) 



a 2^A^2(j_p^ 
271 °^ mz2 iV(l-p2z2) 



-Pz(l-z') 



+ 



4p(l-z2)(4p-z-2pz2) 



5(l + p2+2pz) 



(^Ao(mf^) -Ao(mz^)) 



1 



2mz'z(5 - 4P^ - (3 - 2p2)z2) _ ^p(i _ (2 + pz) 



5 



+ 



5^ 



+ 



RM^R ^ 2mz2(6p3 + (5-8p2+4p4)z + (5-12p2)Pz2 
sP(l + Pz) L 

(3-4PV2PV-(3-4P^)Pz^) 
2p5(l-p2)(2-p2 + pz-z2_p23 

- [8mz'^ + 4P5(P + z)mz^ + s^{2- 2p2 + p^ + 2pz + p¥) 

2mz4z(5-8p2-(3-2p2)z2) 

2mz25P(l + p2 - 2(1 - 2p2)Pz - (1 - p2)z2 + (1 - p2)Pz3) 
^¥(P + z)(l + Pz)lcf 



1 



+ 



1 



[4mz^P(P + z) (l - 3p2 + (1 - 2p2)pz + 2z2p2 + p^, 



5(1 + PZ)(1 + P2 + 2PZ) 

+ 25(1 - 4p2 + p6 + (2 - 10p2 + 3p4)pz + (3 - 5p2)p2z2 
+ 2(3-p2)pV+4pV + pV)]5o(z) 
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- i[4(2mzVmz25p(P + z))(l + pz) 



2 - 2p2 + p4 + 3Pz + 3p V + ^h^) ] (z) 



+ 2sl 



I6mz^ + l6mz\^^{^ + z) 
+ 2mz^/ (2 - p2 + + 2( 1 + 2p2)Pz + (2 + P^)^ V) 



+ 



+ 



85p 



^(l^p2 + 2Pz) (^o^'"^') -^o^'"^'^ 

(16mft2(^^2 _^^2^) _4^^y2^(i _ p2))2(5 _4p2 _ (3 _2p2)22 

+ Smb^s^{2 + (5 - 4p2)pz - 2z^ - (3 - 2p2)pz3) 

+ 52(l_p2)(4p_(5_3p2_4p4)^_4P^2^(3_3p2_2p4)^3)^ 

- g^p^Y^pp^ [(^16m^ (wz, -mw ) -4mH' 5(1 -p )j 
(6p3 + (5 - 8p2 + 4p4)z + (5 - 12p2) pz2 

- (3-4p2 + 2pV-(3-4p2)pz4) 
+ 16mfe2^P(l - p2) (2 - 3p2 + (3 - 2p2)Pz - (1 - 3p2)z2 

- (2 - p2)pz3 - ^h"") + / (I - p2)2 (^8p - 2p3 - (5 - 4p2 - 4p4)z 

- 9Pz2 + (3-4p2-2pV + 3pz^)]< 



+ 



1 



8^ 



(B.13) 



64mt^ + I6mh\s{ 1 - P^) - Smw^ - 2zs^) 
+ 16mw\s{l - p2) +4mb^)+32mb^mw^s{l + pz) 

+ 4sW (3 - 2p2 + p4 + 4(2 - p2)Pz + 2P V) + 8/mw^P( 1 - P^) (P + z) 
+ 5^(1 _ p2) (3 _ 2p2 + p4 + 2(1 + p2)Pz + 2PV)" 
1 

32^ 

z(5-8p2-(3-2p2)z2 



+ 



- 16 (4mb^ + mw'^{s {I - P^) +4mb^) - ^mb'^mw^^ 



\6mb^s (4( 1 + p2)p - (5 - 17p2 + 8p4)z 
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- 4(l-p2)p^2^(3-9p2 + 2pV) 

+ 64mwW*p(l + p2 + 2pz-(l-|3V) 

+ V (8(2 - p2)p3 + (5 _ + 2ip4 - 8fy 



+ 8(2-p2)p3^2_(3_p4_2p6)^3 

+ 8mw V(l - p2) (2(1 + p2)p _ (5 _ p2 _ 8p4)^ 

- 2(l-p2)pz2+(3_p2_2p4)z3) 

+ 5^(1 - (4(1 + p4)p - (5 - 22p2 + 9p4 - 8|36)z 



- 4(1 - 2p^ - p*)pz^ + (3 - 4P^ + 3P* - 2p6)z3 

1 

45(l + pz)(l + p2 + 2Pz) 

P(P + Z) (1 - 3p2 + ( 1 - 2p2) pz + 2P V + p3^3^ 

- Smb^s (1 - 6p2 + 3p6 + (2 - 16p2 + Vp^ + 2p6)pz 

+ (4 - + 3p4)P V + (9 - 4p2 - p4)p3z3 + 2(3 - p2)pV ^ p5^5^ 



C, 



O T O O O \ 

I6mh {mh —mw )—4mw s{l — ^ )1 



+ 



^2(1 + p2 +2pz) (1 - p2) (2 - 5p2 + p4 + (1 _ 2p2 _ 2p4)P^ 
(3 - 2p2) P V + (2 + p2) p3,3 + p4^4^ j 



1 r 

85 



64m^6( 1 + pz) - 16m/ (8m,y2(l + pz) - ^(1 - - (3 - p2)p2 - 2pV)) 



+ mb^i32{2mw'^ + mw^s{l + ^z)){l + ^z) 



2 ^32(2mv7'^ + mw'^i 
+ 4^2(3 - 4pV + 9pz - p'z + 2p3z3 + lOp^z^ - 2p3z - 2p2 + p^; 
+ 1 6mw^s{\ - p2) ( 1 + pz) +Smw^s^^{l-^^){^ + z){l + pz) 

1 



+ ^3 ( 1 _ p2 ) ( 1 + p^) (^3 _ 2p2 + p4 + 2 ( 1 + p2) Pz + 2P V) " 

,56m/ - 12Smb^{6mw^ + + pz)) 

- 64mw^5(l - p2) + 32m^4 (24mw'^ + 2mw^s{3 + p^ + 8pz + 2pV) 

- /(-3p2 + p4 + 1 + 2p4z2 - 5P V - 6Pz + 2p3z)) 

- 16mivV(l-p2)(l+3p2 + 4pz) 

- 8m/ (32mw^ + 16miv^5(l + P^ + 2pz) + 2mw^/(5 + p^ + 12pz + 2(2 + p2)p V) 



37 



dz 



- Am^^s^ ( 1 _ p2) (3 + 3[34 ^ 4( J ^ 2p2) pz + 6p¥) 

- 5^(1 _ p2) (1 ^ p2 ^ 2Pz) (3 - 2p2 + p4 + 2(1 + p2)Pz + 2p V)] (z) 

+ (z^-z), 



47C'"°" mn/2 iV(l-pV 
4P(l-z2)(4p-z-2pz2) / 



1 



+ 



.(l + p2 + 2Pz) v^o^'^^'^-^o^'^^'O 
2miy22(5-4p2-(3-2p2)z2) 

sP (2 - (7 - 8p2)Pz - 2z2 + (3 - 4p2)Pz3^] 5o 
^P^[2.,2(_6p3_(5_8p2 + 4P^), 

- (5 - 12p2)pz2 + (3 _ 4p2 + 2p4)z3 + (3 _ 4p2)pz4 j 
+ 25P(l-p2)(^2 + 3p2-(3-4p2)Pz 

- (l+6p2)z2 + (l_2p2)Pz3 + 2p2z4)]5S' 

- -[8mij^-4mi?25(2-3p2-pz) 
s L 

+ ^2(6- iop2 + 5p4 _ 2(1 _ 2p2)p^ + p222)] 



^p 



2m7y4z(5-8p2-(3-2p2)z2) 

- iMH^S^ (1 + P' - 2(5 - 6p2)pz - ( 1 - p2)^2 + 3(1 _ [32)(3^3^ 
+ 52[32j^p_2|33_(7_i3p2^8p4)^ 

(l-2p2)Pz2_2(l-p2)p2^3^jc^^ 



+ 
+ 



+ 
+ 



5(1 + PZ)(1 + P2 + 2PZ) 

" - 4m^2p(p ^ ^) J _ 3p2 ^ ( J _ 2p2) ^ 2p2z2 + p3^ 
25 (1 - 3p^ + (2 + 2p2 - 5p4 - 4p6)pz + ( 1 + 3p2 - 6p4)p2r 

2(1 + p2)p5^3 +4p6^4 ^ p5^5^ j^H(^) 

(8m//4 - AmH^s{l - 3p2 - pz)) (1 + pz) 



(B.14) 
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+ s\6- 10p2 + + (3 - + 4p4)Pz - (1 - 4p2)pV + p3^3)j ^h^^^ 

+ J-[-16m7y^ + 16m//^5(l-2p2-pz) 
2s L 

- 2mHV(6- 13p2+ 10p4-6(l -2p2)Pz+ (2 + p2)pV) 

- 53(^2 + p2-6pV4p6 + (5-10p2 + 8p4)pz 

+ (l+2p2)pV + pV)]<(z)| + (z--z). (B.15) 
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